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Abstract 

The objective of this work is to optimally design the thermal components of a system that uses 
carbon dioxide (C0 2 ) from the Martian atmosphere to produce oxygen (0 2 ) for spacecraft 
propulsion and/or life-support. C0 2 is thermally decomposed into carbon monoxide (CO) and 0 2 
followed by the electrochemical separation of 0 2 . The design of the overall system and its various 
individual components depends on, among other things, the fraction of the stoichiometric yield 
of 0 2 that can be realised in the system and the temperature of operation of the electrochemical 
separation membrane. The analysis performed indicates that a substantial reduction could be 
obtained in the mass and power requirements of the system if the unreacted C0 2 were to be 
recycled. The report also discusses the concepts of an optimum temperature of the Zirconia cell 
and impracticality of plant operation at low cell efficiencies. The design of the thermal equipment 
would be such that the mass and power requirements of the individual components and of the 
overall system would be optimised. 

Introduction 

There has been tremendous interest in recent years on the need for in-situ-propellant-processing 
(ISPP) [1,2]. More than 80% of a spacecraft’s mass is due to the propellant. Hence production 
of propellants at locations remote from the Earth is essential for frequent and extended space 
explorations. Here, the production of Oxygen from the predominantly C0 2 atmosphere of Mars 
is explored. The system that would be used on Mars would have a filter at the C0 2 intake end to 
remove the dust and particulates from the supply gas. It would also have a liquefaction and 
storage unit at the 0 2 production end. These components are not being designed by UA/NASA 



2 4 3 m 

/<!/&? o 


IB-10 


SERC presently. However, they will be part of our future investigations. The system currently being 
designed incorporates the sub-systems necessary for supplying the C0 2 at Martian pressures and 
temperatures. This system will be designed, fabricated and tested extensively on our ground- 
based facilities. To distinguish it from the future flight-tested system for Mars, the ground-based 
testing system will be referred to as the Test Bed". 

Description of the Test Bed 

The Test Bed, shown in figure 1, can roughly be divided into 4 sub-systems: 

1. The simulation sub-system consists of C0 2 supply and a Cryo-Vacuum Chamber cooled by a 
Cryo-cooler. This simulates the Martian ambient conditions of 6.4mbar pressure and a 
temperature of 200K. Though the Martian atmosphere contains only 95.3% C0 2 , presently the 
supply gas is composed entirely of C0 2 . 

2. The compressor, heat exchanger and heater form the C0 2 preparatory sub-system. In order 
to avoid a very bulky system at high C0 2 mass flow rates, the supply gas is pressurized from 6.4 
mbar to higher pressures. At the present time C0 2 is assumed to be pressurized to ibar. The 
reasons for this are two-fold. The Test Bed is operated on the Earth and operating the system at 
1 bar minimises the probability of gas leaks. Secondly, electrochemical separation membranes 
have been tested extensively with other gases at 1 bar and higher pressures and its performance 
at lower pressures is not known currently. The effect of lower pressures on the electrochemical 
separation is the subject of a related investigation at UA/NASA SERC. After the compressor, the 
pressurized C0 2 enters a waste heat recovery heat exchanger to gain energy from the gases 
exhausting from a Zirconia (Zr0 2 ) Cell, which will be discussed later. The C0 2 supply gas then 
enters a heater where it reaches the same temperature as the Zr0 2 electrochemical separation 
membrane unit. Thermal decomposition of C0 2 occurs in the heater. Since the temperature of the 
gases entering the separation unit is the same as that of the unit itself, the separation unit is not 
subjected to any thermal shocks. 

3. The mixture of CO, 0 2 and C0 2 passes through the Zr0 2 cell, which electrochemically 
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separates the 0 2 from the mixture. This 0 2 is then fed to a storage device for later use. The cell 
forms the oxygen separation sub-system. 

4. The exhaust from the Cell, after passing through the heat exchanger, loses energy in a radiator 
before passing through a polymeric membrane separator which separates the unreacted C0 2 . 
The radiator is essential because polymeric membrane separators cannot withstand the high 
exhaust temperatures. There is a substantial pressure drop across the membrane separator, 
which requires a second compressor to re-pressurize the C0 2 to 1 bar. The C0 2 then re-enters 
the loop just before the heat exchanger. The system is depicted in figure 1 . 

Figure 2 illustrates the system without recirculation. The exhaust gas, after passing through the 
heat exchanger, is expelled from the system. Henceforth, the system with recycling is referred to 
as case 1, and the one without, as case 2. The number of moles of CO, C0 2 and 0 2 flowing at 
various points in the system for both cases 1 and 2 is shown in table 1 . 

Nomenclature 

e = Fraction of 0 2 produced that is electrochemically separated by the cell. 

(i = Membrane Separation Factor = Amount of C0 2 at Point 9 
sep Amount of C0 2 at Point 8 

This is only for case 1 . 

a = Cell Factor = moles of unreacted C0 2 for each mole of 0 2 produced 

This determines the cell efficiency i cell as: 

t cell = moles of Q-. actually produced = 2/(2+ a) 

” e theoretical max. no. of moles of 0 2 possible 

The Cell reaction is: 

[(2+a)/e]C0 2 = [2/e]CO + [1/e]0 2 + [a/e]C0 2 
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Fig 2: Test Bed Without Recycling of Carbon Dioxide 
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Table 1 : For Oxygen Production of 1 mole 


Point 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 


Case 1 

(2+a.-a sep tt)C0 2 

G 

(2±o)C0 2 

e 

(2+a)C0 2 

€ 

(2±ja)C0 2 

€ 

(2+a)C0 2 

e 

[2/g]CO+ [(1 -€)/e]0 2 + [a/ € ]C0 2 

[2/€]C0+[(1-€)/€]0 2 +[a/€]C0 2 

[2/€]C0+[(1-€)/€]0 2 + [a/e]C0 2 

fPsep a/€ l C0 2 

[Psep a/e l C0 2 


Case 2 

(2±o)C0 2 

€ 

(2+«) CQ 2 

g 

(2±a)C0 2 

€ 

(2±a)C0 2 

€ 

(2+a)C0 2 

€ 

[2/€]C0+[(1«€)/€]0 2 + [a/e]C0 2 
[2/e]C0+[(1-e)/e]0 2 + [a/G]C0 2 


Assumptions 

The thermodynamic analysis has been performed based on the following assumptions: 

1. The Membrane Separator filters out all the CO and 0 2 from the exhaust gases. It is known that 
the purity of the separated gases is very high. Hence this is a valid first-order approximation. 

2. Only a steady state analysis has been performed. 

3. Pressure drop and heat losses in the pipes, valves and bends are negligible. 

4. The Compressors have been assumed to be isentropic, with specified mechanical and 
isentropic efficiencies. Radiation losses from their surfaces have been neglected since they are 
small in comparison. 

The analysis gives the total heat transferred from the exhaust to the fresh C0 2 . Using an iterative 
procedure, the cold side outlet temperature and the Logarithmic Mean Temperature Difference 
(LMTD) are determined. Knowing the mass flow rates, the heat exchanger can be designed. The 
power requirement of the Zr0 2 cell - the Nernst potential and the ionic component of the current 
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corresponding to the oxygen ion flow can be calculated precisely. For a fixed 0 2 production rate, 
it is a constant. 

Discussion of Results 

The analysis was performed with the help of an interactive program written in Pascal. The 
program has pull-down menus using which the user can change any of the input variables. This 
was necessary because the optimal cell temperature, cell efficiency and membrane separation 
factor is presently not known. Concurrent research is being done at UA/NASA SERC to investigate 
the optimum system performance conditions. Some graphs have been obtained from the data 
generated by using the program for an oxygen production rate of 10 kg/day, which corresponds 
to the production rate needed for an unmanned sample return mission [2]. Figure 3 shows the 
variation of the specific heats of C0 2 , CO and 0 
have been used to determine these curves. 

This variation in the specific heat (Cp) precludes 
treating C p as a constant, say, across the heat 
exchanger. Over the entire range of 
temperatures considered, it is seen that C p 
increases with temperature for all three gases, 
with the curve for C0 2 increasing most sharply. 

Figure 4 shows the graph of heater power with 
cell temperature (T 5 ). The power consumed by 
the heater increases with T 5 as expected, for 
low values of Te. It reaches a maximum, and then, contrary to expectations, starts dropping. 
Figures 6 and 7 show plots of temperature of the cell (T 5 ) with the outlet temperature T 4 of the 
cold C0 2 from the heat exchanger. Though not apparent from the figures, (T 5 -T 4 ) actually 
increases with T 5 . The anomalous behaviour of figure 4 can be explained by the fact that though 
C p itself increases, and fT 5 -T 4 ) also increases, the integral of C p dT over the temperature range 
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Fig 4: Heater Power - Case 1 


Fig 5: Heater Power - Case 2 




Fig 6: T 4 v/s T 5 - Case 1 Fig 7: T 4 v/s T 5 - Case 2 


T 4 to T 5 actually decreases after a point. This is opposed to the almost linear increase of heater 
power with temperature in case 2 as shown in figure 5. In this case, though the mass flow rates 
of C0 2 at points 4 and 5 are the same as in case 1 , T 3 , and hence T 4 , are lower for the same 
value of T s . This is because the recirculating C0 2 is at a higher temperature than that of the C0 2 
at point 2. Hence, the increased (T 5 -T 4 ) term counters the decrease in the slope of the C p curve 
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enough to cause the heater work to increase. Here, and in discussions later in this report, we 
freely use the '(T 5 -T 4 )' argument. Although the empirical relations used for C p were fourth-order 
polynomials in temperature, terms such as (T 5 *-T 4 4 ) etc. depend essentially on (T 5 -T 4 ). Note that 
the magnitude of heater work is greater for case 2 than for case 1, for a fixed T 5 . It may be 
expected that at cell temperatures higher than those considered, figure 5 behave similar to figure 
4. However, in practice, cell temperatures higher than 1 500K are not of interest. 




iou.u - 





-299.04 



-279.04 

140.0 - 




\ 

- 259.04 

o 

■o 

\ 

L 239.04 E 

\ 1 20.0 - 

\ 


CD 

\ 

i-21 9.04 


\ 



\ 

-199.04 -2 

9 100.0 - 

\ 

O 


\ 

-179.04 02 

— 80.0 - 

\ 

-159.04 | 

o 

\ 

-139.04 U " 

O 

60.0 - 

\. 

-119.04 O 


N. 

-99.04 

• 



it 40.0 - 


-79.04 



-59.04 



- 39.04 

20.0 - 

1 1 1 1 1 i r 


0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0 . 

8 


Cell Efficiency 




Fig 8: C0 2 Intake - Case 1 Fig 9: C0 2 Intake - Case 2 


Figures 8 and 9 show the C0 2 intake required for cases 1 and 2 respectively. Clearly, the intake 
required increases sharply as the cell efficiency approaches zero. The intake in cubic-feet per 
minute (cfm) at the inlet pressure of 6.4 mbar gives us an idea of the size of pumps that we will 
require. Note that for a given cell efficiency, the intake required in case 2 is greater than that in 
case 1, especially at lower efficiencies. Since the slope of the curve decreases rapidly with 
increasing efficiency, and keeping in mind that the intake determines the capacity of the 
Compressor C v the conclusion is that the compressor capacity becomes impractical at very low 
efficiencies. For a start, we take an efficiency of 0.2 which is reasonable at this point of time. The 
pumping speed required is about 280 cfm for case 2 and about 115 cfm for case 1 . If we take an 









IB-18 


off-the-shelf vacuum pump which can fulfil these requirements, then for case 2, an S Series Rotary 
Vane Pump Model S630F manufactured by Leybold Vacuum Products Inc. will suffice. It weighs 
1450 lb and has a power requirement of 25 hp. For case 1, where the C0 2 is recycled, a much 
smaller pump such as the Leybold Model S160C, which weighs 310 lb and consumes 7.5 hp can 
do the job. Also required in case 1 is a much smaller pump (compressor C 2 ) to compensate for 
the pressure drop across the membrane separator. For this a small pump like the Leybold TRIVAC 
Rotary Vane Pump Model D2A will do. This pump weighs 41 lb and consumes 0.33 hp. This 
clearly indicates that at this point, though case 2 appears compelling from the points-of-view of 
system simplicity and reliability, we must recycle the C0 2 in order to design a reasonably compact 
system. Note that we must also consider the mass added to the system by the radiator and the 
membrane separator. However, the reduction in compressor mass and power obtained by 
recycling offsets the increase in mass due to the radiator and membrane separator. Having 
decided on recycling, we take an optimistic look at what our requirements would be if we were 
successful in obtaining a cell efficiency of, say, 0.5, which gives us a CO z intake of about 70 cfm. 
To take care of this requirement, a first approximation would be a pump like the Leybold S100C, 
with a mass of 220 lb and a power requirement of 5.0 hp. This makes the system design very 
feasible as far as its mass requirements go. Obviously, for an actual mission, a much more 




Fig 10: Heater Power v/s t cell 
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sophisticated pump with better materials and design would be used. But it would depend on 
whether an efficient Zirconia cell can be designed. Figure 10 is a graph of heater power versus 
cell efficiency for case 1 . In the light of information obtained from figures 4 and 1 0, we realise that 
at low efficiencies the heater power increases because of increased C0 2 mass flow rate. But the 
interesting point made by this graph is that at higher cell temperatures and lower efficiencies, the 
effect of the cell temperature prevails, and causes the slope of the curve to decrease.ln fact, the 
curve for a cell temperature of 1400K is almost flat, suggesting that there could be a temperature 
T 5 between 1400K and 1200K at which the load on the heater is almost constant. Also important 
is the fact that it is a lower temperature difference (T 5 -T 4 ) that "overcomes" the increase in the 
slope. Figure 1 1 , for case 2, is clearly a follow-up of figure 5. Here, the temperature difference (T s - 
TJ does not fall sufficiently to cause the slope to decrease. 

Plans For The Future 

1. The system design depends, to a large extent, upon the system operating pressure. The 
optimum system pressure needs to be examined. If the system pressure were low, then the option 
of using a blower (lower mass) instead of a compressor will be examined. 

2. The heat exchanger design presently limits the exhaust gas temperature to 700K, due to fears 
of carbon deposition. Suitable changes will be made in the design when a related investigation 
by Prof. D.C. Lynch, Department of Materials Science, University of Arizona, provides results on 
the seriousness of carbon deposition. 

3. The pressure and heat losses in the system components and piping will be incorporated in the 
analysis and 

4. The individual components will be designed and tested. 

Acknowledgements 

The graphics routines of the Pascal program used in the analysis were written by Mario Rascon 
and Kirby Hnat. Their help is gratefully acknowledged. 


IB-20 


References 

1. K. Ramohalli, E. Lawton & R. Ash, "Recent Concepts in Missions to Mars: Extraterrestrial 
Processes", Journal of Propulsion and Power, Vol 5, pp 181-187, 1989 

2. R.L Ash, W.L Dowler & G. Varsi, ’ Feasibility of Rocket Propellant Production on Mars", Acta 
Astronautica, Vol 5, pp 705-724, 1 978 

3. Thermophysical Properties of Matter, The TPRC Data Series, Vol 6, p 50, 145 & 1 54, IFI/Plenum, 


1970 


